Miniature plasmonic wave plates 



O 

(N 

IX, 
m 



CZ3 

O 

-f— » 

q 

O 
Oh. 



> 

m 
o 

O 

o 



X 



Aurelien Drezet, Cyriaque Genet, and Thomas W. Ebbesen 
Isis, Louis Pasteur University, 8 allee Gaspard Monge, 67000, Strasbourg, France 

(Dated: February 3, 2010) 

Linear birefringence, as implemented in wave plates, is a natural way to control the state of 
polarization of light. We report on a general method for designing miniature planar wave plates 
using surface plasmons. The resonant optical device considered here is a single circular aperture 
surrounded by an elliptical antenna grating. The difference in short and long axis of each ellipses 
introduces a phase shift on the surface plasmons which enables the realization of a quarter wave 
plate. Furthermore, the experimental results and the theoretical analysis show that the general 
procedure used does not influence the optical coherence of the polarization state and allows us 
to explore completely the surface of the unit Poincare sphere by changing only the shape of the 
elliptical grating. 

PACS numbers: 42.25.Lc, 42.70.-a, 73.20.Mf 



Surface plasmon polaritons (SPPs), electromagnetic sur- 
face waves existing at the interface between a dielec- 
tric and a metal are particularly sensitive to tiny 
variations in their local electronic environments. This 
creates new opportunities and applications for photon- 
ics Q by simply texturing a metal surface. For example, 
metal films structured with two dimensional subwave- 
length hole arrays present remarkable properties such as 
the extraordinary optical transmission (EOT) which is a 
clear signature of SPP-light interaction [3|45| . In this par- 
ticular context, several studies have started to address 
polarization issues, discussing in this respect the influ- 
ence of the individual hole shapes. Elliptical or rectan- 
gular apertures can behave like polarizers, following the 
Malus law of absorption (see also ref. Q for similar 
work on elliptical nanoparticles) . However these struc- 
tures do not display linear birefringence. Linear birefrin- 
gence is absolutely central in optics since it allows full 
control of the state of polarization (SOP) of light with- 
out absorption. A half wave-plate rotates the plane of 
polarization while a quarter wave-plate converts linear 
polarized light into a circular one, the combination of 
the two enabling a complete exploration of all polariza- 
tion states. 

In this letter, we report for the first time, both ex- 
perimentally and theoretically, the design and charac- 
terization of a plasmonic optical wave-plate. In order 
to obtain the linear birefringence, we have developped a 
modified version of the single circular nano-aperture sur- 
rounded by periodic circular corrugations, also known as 
a bull's eye structure (lflj. Such an optical grating acts 
as a miniature antenna presenting huge EOT for optical 
wavelength inside a narrow band centered on the SPP 
resonance 1^, U | ■ The specificity of the structure pre- 
sented here is its unique ability to control the SOP of the 
electromagnetioc field going through the aperture. This 
is achieved by introducing a well defined excentricity in 
the grating geometry which in turn modifies the phase 
of the excited SPP and consequently the polarization of 



the transmitted light. This resembles in a wide sense the 
phase matching in distributed feedback lasers (DFB). To 
fully characterize the optical behaviour of our device, a 
genuine polarization tomography of the isolated subwave- 
length aperture had to be implemented. Furthermore we 
have developped a microscopic (dipolar) model to link 
structural design with change of SOP. 
Beside its fundamental interest, such control over the 
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FIG. 1: (a) Scanning electron microscopy image of the el- 
liptical bull's eye structure. The scale bar is 2/xm long. The 
green arrow is an excited point dipole Pm located in M. SPP 
(red arrow) launched in M propagates along the radial direc- 
tion MO to reach the hole located in O. w is the angle 
between MO and Pm- (b) Comparison between the exper- 
imental white light transmission spectrum (red curve) asso- 
ciated with the structure shown in (a) and the theoretical 
prediction (black curve) obtained with the 2D dipoles model. 
The transmission scale is dimensionless and corresponds to 
the ratio between the transmitted and the incoming power at 
the level of the aperture. 

SOP can be used broadly in photonic applications requir- 
ing local addressing, e. g., detectors [12j, displays [ll| and 
compact circular polarization antennas [l3| . In addition, 
ultrafast opto-magnetic data storage has been demon- 
strated with femtosecond lasers and circularly polarized 
14). The device demonstrated here creates both 
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the right helicity and the large fields in a tiny volume 
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favorable for such purposes. 

For our experiments we consider a bull's eye structure 
made of 8 grooves and fabricated by FIB milling in a 
300 nm thick Au film (Fig. 1). The hole diameter is 
260 nm and the grooves width and depth are 370 and 
80 nm respectively. The groove shape is chosen to be 
elliptical with the long axis a n = n ■ P + P/4 and the 
short axis b n = n ■ P. Here P — 760 nm is the pe- 
riod of the grating (which equals the SPP wavelength 
Xspp for a laser excitation at 785 nm [lH) and n is an 
integer going from 1 to 8 (see Fig. 1(a)). Also shown 
in Fig. 1 (b) is the transmission spectra of the structure 
with a resonant peak at Ao — lllnm. The measured 
extraordinary transmission efficiency (larger than one) is 
a direct signature of the involvement of SPP [2]. The 
presence of this transmission peak proves that despite 
the small increment of SL = P/4 between the long and 
short axis of the ellipses the structure still behaves like 
a miniature antenna. We can justify our choice for the 
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FIG. 2: (a) Sketch of the optical polarization tomography 
setup. P, A/4, A/2 are respectively polarizers, quarter wave 
plates, and half wave plates located in the input and output 
beams. The image is recorded with a CMOS camera and the 
light source S is a laser diode emitting at A = 785nm. L\ 
and Z/2 are two objective lenses (x50, NA = 0.55) and (x40, 
NA — 0.6) respectively, (b) Typical camera image of the light 
transmitted through our structure (the scale bar is 2 (im) . (c) 
Cross cut of the intensity profile along the white dotted line 
shown in (b). 

grating symmetry on theoretical grounds. In our model 
we discretize the grooves into a sum of point dipoles P m 
proportional to the local electric field at M. Each dipole 
is excited coherently by the light impinging normal to 
the metal film and SPPs are launched in the direction of 
the central nanohole where they excite an in-plane radi- 
ating dipole [16(. To reproduce completely the system, 
we introduce a second transmission channel in which the 
central dipole is excited directly by the incident light. 
The interference between these two channels leads to a 



Fano like effect 17J resulting in the observed transmis- 
sion peak. The relative (complex) amplitude between 
these two channels was fitted to reproduce the spectra of 
Fig. 1(b) [18| . The good agreement between our model 
and the data (see Fig. 1(b)) allows us to use it for pre- 
dicting the optical behavior of the structure at a given 
A. The principle of the device can be illustrated by con- 
sidering only the point dipoles located along the short 
and long axes of the ellipses. It is thus clear that SL 
corresponds to a phase shift 4>spp — ^ttSL/Xspp = n/2 
between SPPs propagating along the long (y) and the 
short (x) axes. Additionally the coupling between the 
incident light and SPPs depends on the cosine of the an- 
gle w between the radial vector MO and P m (Fig. 1 
(a)). It means that if the incident linear polarization is 
switched from a direction parallel to the x axis to a di- 
rection parallel to the y axis then the radiating central 
dipole will change from ax (where a is a constant) to 
e !X / 2 tty. From the point of view of this idealized picture 
(which neglects damping) we deduce that the system be- 
haves like a birefringent biaxial medium, i. e., a perfect 
quarter wave plate, with fast and slow axes parallel re- 
spectively to the x and y axes. Obviously if we now take 
into account all the dipoles as well as the Fano interfer- 
ence effect and the finite value of the SPP propagation 
length Lspp (damping) in the structure the actual result 
will naturally deviate from this idealized case 

In order to study experimentally the SOP conversion 
by our structure, we carried out a complete polariza- 
tion tomography (2]J using the optical setup sketched 
in Fig. 2(a). A laser beam at A = 785nm is focused 
normally onto the structure by using an objective L\. 
The transmitted light is collected by a second objective 
Li forming an Airy spot on the camera (see Fig. (2b-c)) 
as expected since the hole behaves like a point source 
in a opaque gold film. In our experiments, the inten- 
sity is thus defined by taking the maximum of the Airy 
spot shown on Fig. 2(b). The SOP of light is pre- 
pared and analyzed with half wave plates, quarter wave 
plates, and polarizers located before and after the ob- 
jectives [21.-23]. We recall that the complete knowledge 
of the SOP require 6 intensity projection measurements 
(ijfj) = (|E-ai| 2 ) made along the 4 linear polarization vec- 
tors x, y, p = (x + y)/y/2, m = (x - y)/ \/2, and along 
the two circular polarization vectors L = (x + iy)/v2), 
£y)/v2)- It is convenient 22 1 to introduce the 
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goal of this polarization tomography is then the deter- 
mination of the 4x4 Mueller matrix M. characterizing 
the transformation of the input Stokes parameters dur- 
ing the interaction of the laser light with the structure. 
In order to write down the full Mueller matrix, we mea- 
sured 6x6 intensity projections corresponding to the 6 
previously mentioned unit vectors for the input and the 
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output polarizations 24 j. 

At first, the isotropy of the bare setup was checked 
by measuring the Mueller matrix v\/^ ass with a glass 
substrate. Up to a normalization constant, we deduced 
that _A/f§^ ass is practically identical to the identity ma- 
trix X with individuals elements deviating from it by no 
more than 0.02. It implies that the optical setup does 
not induce depolarization and that consequently we can 
rely on our measurement procedure for obtaining M.. 
Optical depolarization (i. e, losses in polarization coher- 
ence) can be precisely quantified through the degree of 
purity of the Mueller matrix defined by 25] D (M) = 

( Tv[M l%l Ml ° ) V2 < 1. We find D (X§ lass ) = 0.9851. 
We impute the residual depolarization (1 — D ~2%) to 
the lenses and to alignment errors. It should be noted 
that the incident illumination spot size on the sample 
was varied between 2 and 20 pm without affecting the 
matrix, i. e., without introducing additional depolariza- 
tion. In the following experiment, we consider the case 
of a large gaussian spot with FWHM=20 /im in order to 
illuminate the whole structure [2q |. We then measured 
the Mueller matrix of our structure and found that: 



^exp. = 



/ 1.000 
0.111 



0.107 



0.972 

0.004 -0.002 
\ 0.001 -0.017 



-0.008 0.000 \ 
-0.002 -0.004 
0.306 -0.932 
0.934 0.294 / 



(1) 



which is clearly block diagonal, up to experimental errors. 
It is also remarkable that we have D(A^ exp- ) = 0.981. 
This means that despite the existence of the SPP trans- 
mission channel, the polarization coherence is not lost 
during the propagation through the structure. This situ- 
ation contrasts with previous SOP tomography measure- 
ments on metallic hole arrays in which the polarization 
degrees of freedom were mixed with spatial information 
responsible for SPP-induced depolarization (StJ. Beside 
these two points, the matrix A4 ex P" exhibits several in- 
teresting symmetrical features which relate to the po- 
larization properties of the device. First, it can be ob- 
served that in our experimental procedure the polariza- 
tion in the Airy spot (see Fig. 2(b)) is homogeneous [28j . 
This means that in our analysis we are actually doing the 
polarization tomography of the central radiating dipolc, 
i. e., we are dealing only with the SU(2) point symmetry 
of the Mueller matrix. In such context, the rectangular 
point symmetry group C 1 ?,, of the ellipse imposes that 
the (2 x 2) Jones matrix [22| J connecting the incident 

electric field \ E*P-, E}J^] to the transmitted electric field 
must be diagonal in the x and y basis, i. e 

3 ( n o I where /3 = pe 1 ^ is a complex number 
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analogy with bulk optics, p and <p measure respectively 
the relative dichroism (i. e. the relative absorbtion) and 
the birefringence of this biaxial 2D medium. Clearly </> is 



reminiscent of <pspp discussed above. Using J we obtain 
the theoretical Mueller matrix 
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2Re {/3) 
21m 1/3) 






-21m (/3) 
2Re (/?) 



(2) 



which is similar to .M ex P- and in particular satisfies the 
symmetries Moi = Mio, Moo = Mn, M22 = M33, 
and M23 = — M32 observed experimentally. We de- 

1 /2 

duce P = ( mZ+MoI ) and tan = M 3 2/M 2 2. Using 
Eqs. 1,2 we obtain p ~ 0.898 and 4> ~ 72.5°. Recipro- 
cally by injecting the previous values for p and cf> in M 2 " 
the result do not differ from .M ex P- by more than 2%, 
in agreement with the value obtain for the residual de- 
polarization. Then, using the fitting parameters already 
considered in the transmission spectrum Fig. 1(b), we nu- 
merically calculate the Mueller matrix predicted by the 
2D dipole model and obtain 



M 



2D 



( 1.000 0.089 0.000 0.000 ^ 

0.089 1.000 0.000 0.000 

0.000 0.000 0.446 -0.890 

\ 0.0000 0.000 0.890 0.446 / 



(3) 



which is close to M C2v and M ex P- and corresponds to 
D (M 2D ) = 1. This numerical model is sensitive to small 
variations of the fitting parameters and the agreement 
with the experiment could be probably improved by go- 
ing beyond the paraxial approximation for the incident 
light 0. 

Finally, we considered more closely the consequence 
of the transformation defined by .A/( ex P' by varying the 
linear polarization 9 of the input state every 10° from 
-90° to +90°. On Fig. 3(a) we show the transmitted in- 
tensity analyzed along the 6 fundamental polarizations 
x, y, p, m, L and R as a function of 9. The interfer- 
ence fringes observed are compared with the predictions 
given by the 2D dipole model (dotted curves) and with 
the intensity deduced from the Mueller matrix .M ex P- 
(continuous curves) . In both cases the agreement is very 
good showing once again the consistency of the different 
measurements and deductions. Furthermore, this can be 
geometrically illustrated by using the Stokes vector 22 1 
defined by S = [Six + S^y + Sazj/So- The surface drawn 
by the input Stokes vector is called a Poincare sphere and 
has the radius D = 1. As shown in Fig. 3(b) the operator 
_A/f ex P- defines a geometrical transformation connecting 
this Poincare sphere to an output surface with a char- 
acteristic radius D (A4 ex P'). This experimental surface 
is very close to the ideal sphere D = 1 in agreement 
with the absence of net depolarization as discussed ear- 
lier. The experiment shown on Fig. 3(a) is also repre- 
sented on this sphere. From Eq. 2, we deduce that if 
the input Stokes vector explores the equator correspond- 
ing to linear polarizations, then the output Stokes vector 



4 




angle (degree) 

FIG. 3: (a) SOP analysis of the output beam for a linearly 
polarized input beam. The polarization angle is measured rel- 
atively to the x axis. Data points are compared to Eq. 1 (con- 
tinuous curves) and to predictions of Eq. 3 (dotted curves). 
The colors black, green, orange, cyan, magenta, red and blue 
correspond respectively to ^otal' Iz, I y , Ip, I&, Jf,, (b) 
Image of the input Poincare sphere through the transforma- 
tion yV( ex P' (yellow sphere). The black and red data points 
are respectively the input and output Stokes vectors associ- 
ated with the experiment shown on (a). The circle represents 
the predictions deduced from M ex P' . The device also con- 
verts a input L state (black) star into a state shown by a red 
star. 

draws a circle of radius D ~ 1 which is contained in the 
plane y/z = M33/M32 making the angle 90° — 4> = 17.5° 
with the z axis. These predictions are directly consis- 
tent with the observations. Additionally, for a pure left 
circular input SOP, we experimentally obtain S ex P" = 
0.1123x- 0.8984y + 0.3104z in agreement with the value 
deduced from M c ^: S C2 « = 0.1067x-0.9272y+0.2949z. 

To conclude, all this experimental and theoretical anal- 
ysis demonstrate that we have a clear understanding of 
the SPP structure considered here. First, we have p ~ 1 
which implies that the system acts essentially as a bire- 

fringent medium with Jones Matrix J ~ ^ J ^ , i. e., 

a wave plate. Second, the value obtained for <f> shows that 
the system differs slightly from an ideal quarter wave- 
plate for which <fi — 90°. From the point of view of the 
Poincare sphere, this angle measures directly the incli- 
nation of the output circle shown on Fig. 3(b). For a 
perfect quarter wave plate this circle will go through the 
poles, i. e., a complete conversion from linear to circu- 
lar polarization will become possible if the input SOP is 
polarized along p or m. In this context, numerical calcu- 
lations with the 2D dipoles model show that by changing 
slightly the value for the long axis increment SL we can 
change the phase <p continuously. This means than with 
such SPP device we can in principle tailor and generate 
any kind of SOP conversion on the Poincare sphere go- 
ing from the equator (linear polarization) to the poles 
(circular polarization) or vice versa. We expect that the 



SPP control over the polarization presented in this letter 
could have many applications in photonics and in infor- 
mation storage technology. 
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